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Abstract The model of a radial minor merger proposed by lU, which successfully 
reproduces the observed regular shell systems in shell galaxies, is ideal for a test- 
particle simulation. We compare such a simulation with a self-consistent one. They 
agree very well in positions of the first generation of shells but potentially important 
effects - dynamical friction and gradual decay of the dwarf galaxy - are not present 
in the test-particle model, therefore we look for a proper way to include them. 



Tides and dynamical friction in test particle simulations 

We model the luminous and dark matter components of a host giant elliptical (gE) 
and a dwarf elliptical galaxy by analytical potentials. In the simplest model, we 
assume the dwarf galaxy, filled with millions of test particles, to be ripped apart in- 
stantly when it comes close to the center of the gE galaxy. Its stars begin to oscillate 
in the potential of the host galaxy and produce shells at their turning points. 

Such a setup allows us to use large numbers of particles and so to gain sufficient 
contrast to detect all the shells in the simulation, also to investigate the kinematic 
footprint in spectral lines (see Q and 15]) and explore a large parameter space. This 
would be very time consuming for large sets of self-consistent simulations. 

Surprisingly, the agreement with a self-consistent simulation (for more details 
see lHI) turns out to be very good especially in the positions of shells (see Fig. [T] 
and fT\). But the simple model does not involve effects like dynamical friction and 
gradual decay of the dwarf galaxy, so that it cannot simulate phenomena seen in 
self-consistent simulations: the next generation of shells (see H]) and lowering of 
brightness of shells. We thus look for a middle way, where we can still have the 
large contrast available through the use of test particles, yet include some of the 
more complicated effects to make it more realistic. 
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Fig. 1 The graph shows the comparison of histograms of radial distances of shells' particles (cen- 
tered on the host galaxy) in the self-consistent (black) and test-particle (grey) simulations at two 
different time steps. Time equal zero corresponds to the passage of the dwarf galaxy through the 
center of the host galaxy. Both simulations have the same initial conditions. Notice that the bright- 
ness of the outmermost shell (at 280 kpc at 2.55 Gyr) is supressed in the self-consistent simulation. 
This effect has been successfully simulated in the improved test-particle simulations; see |4l- 

In this improved test-particle simulation we use our version of enhanced Chan- 
drasekhar formula with variable Coulomb logarithm to include dynamical friction, 
and we also introduced a gradual decline of the mass parameter of the dwarf galaxy 
potential to better imitate the evolution of tides. For details, see [4 \ and [SJ. 
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